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Abstract 

The twist-2 heavy flavor contributions to the polarized structure function g2ix,Q^) are 
calculated. We show that this part of g2{x,Q'^) is related to the heavy flavor contribution 
to gi{x, Q"^) by the Wandzura-Wilczek relation to all orders in the strong coupling constant. 
Numerical results are presented. 



1 Introduction 



Heavy flavor contributions to polarized and unpolarized deep inelastic structure functions due 
to charm and bottom-quark production cause different scaling violations if compared to those 
due to massless partons. Moreover it is known that in certain phase space regions of x and Q"^ 
these contributions can be large [1]. Reliable determinations of the QCD scale Aqcd therefore 
require a careful account of the heavy flavor contributions. In the case of the unpolarized 
structure functions F2^l{.x,Q'^) and xF^^x.Q"^) the twist-2 contributions with QQ-final states, 
Q = c,b were calculated to next-to-leading order (NLO) [2, 3]. For polarized deep-inelastic 
scattering only the leading order corrections are known for the structure function gi^XjQ"^) at 
present [4]. All corrections mentioned above have been calculated using mass factorization. As 
shown in [5-8] this method fails in the case of the polarized structure functions which emerge for 
transverse polarization already in the massless quark limit. The reason for this lies in a violation 
of covariance due to the omission of contributions oc S.k, where S denotes the spin vector of the 
nucleon and k the parton 4-momentum. 

In Refs. [6, 7] it was shown that the quarkonic contributions to the transverse structure 
functions can be correctly obtained using the covariant parton model. The results in this ap- 
proach are identical to those derived in the local light-cone expansion, cf. [9, 8], for massless 
quarks on the level of the twist-2 contributions. In particular the twist-2 part of g2{x,Q'^) is 
obtained by the Wandzura-Wilczek relation. The effect of quark masses for quarkonic matrix 
elements was further investigated in Ref. [10] using the method of [11]. Also in this case the 
Wandzura-Wilczek relation was found to hold, irrespective of the values of the quark masses 
chosen. In [12, 10] the target mass corrections to the quarkonic contributions to g2{x,Q'^) and 
the other four polarized structure functions [10] were studied. Again the Wandzura-Wilczek 
relation was found to hold relating the twist-2 contributions of gi{x,Q^) and g2{x,Q'^). In a 
more general approach the amplitudes which contribute to the quarkonic matrix elements for 
deeply virtual non-forward scattering were investigated in [13]. In the generalized Bjorken limit 
the non-forward Compton amplitude is expressed by operator matrix elements of vector op- 
erators. However, one usually parameterizes it in terms of scalar operator matrix elements to 
which the (generalized) parton distributions correspond. Therefore relations between the oper- 
ator matrix elements of the vector and scalar operators are implied, which form the origin of 
the Wandzura-Wilczek relations and other integral-relations [8, 10]. The absence of the typical 
integral terms e.g. in the Callan-Gross [14] relation is merely the exception and caused by a 
cancellation of the former, see Ref. [15]. For the polarized case it was shown that the Wandzura- 
Wilczek relation holds for the twist-2 contributions to the respective non-forward amplitudes, 
cf. [15]. This result could be generalized allowing multiple meson production in the final state 
in [16]. Also semi-inclusive processes such as diffrative scattering have been investigated w.r.t. 
the emergence of Callan-Gross and Wandzura-Wilczek relations [17]. Although the variables 
change, the Wandzura-Wilczek relation relates the twist-2 contributions of gf^^'^lXjXpjQ'^) to 
g2^^^ {x , X p , Q"^) independently of the value for xp. Integral relations of similar kind were also 
established for other structure functions and matrix elements [18] and for the fermionic twist-3 
contributions in [10]. 

In the present paper we calculate the heavy fiavor twist-2 contribution to the structure 
function g2{x,Q'^). We show that the Wandzura-Wilczek relation holds also in this case and 
present numerical predictions. 
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2 Longitudinal Gluon Polarization 



The hadronic tensor of the polarized part of the eA'^-scattering cross section in case of pure 
photon exchange^ reads 



A CO- 



p.q 



q^ip.qS" - S.qp") 
{p.qY 



So- denotes the nucleon spin vector, p the nucleon momentum, and q the vector of the 4- 
momentum transfer, with = — and x = Q'^/{2p.q). The polarized part of the scattering 
cross section for longitudinal nucleon polarization Sl, integrated over the azimuthal angle 0, is 



dV(A, ±Sl) 
dxdy 



a 



-2Xy (2-y- ^-^] xg,{x, Q') + SA^^^^l^, Q') 



•(2) 



Correspondingly, for transversely polarized nucleons one obtains 

rfV(A, ±St) 



dxdyd(f) 



±S^2j — 



a 



Q 



xy 



X 



cos(x - 0) 



-2Xyxgi{x, Q^) - AXxg2{x, Q^) 



(3) 



Here M is the nucleon mass, S the cms energy, a the fine structure constant, y = 2p.q/S, A is the 
degree of lepton polarization and St the degree of hadronic transverse polarization, x denotes the 
azimuthal angle associated with Sj^, (10), and gi{x, Q^) and g2{x, Q^) are the polarized structure 
functions which contribute in this case. 

The heavy flavor contributions to the longitudinal and transverse differential scattering cross 
sections are obtained in calculating the corresponding contributions to the polarized structure 
functions gi{x, Q"^) and g^^x^Q"^). 

Let us consider the sub-system hadronic tensor for photon-interactions with on-shell initial 
state partons 



w 



q.k 



'HI/ pa 



(4) 



In assuming that the gluon is longitudinally polarized, i.e. parallel to the proton and parton 
4-momentum 



(5) 



with k — zp, the sub-system hadronic tensor of the polarization asymmetry (4) receives purely 
longitudinal contributions, since the term oc g2^^^°^ vanishes. Clearly (5) is a special model 
assumption which does not describe the general case being discussed in section 3. However, it is 
possible to derive in this approximation the correct expression for the coefficient functions C^*^ 

at twist-2 which contributes to the structure function g^'^^x^Q'^). At leading order in the 

structure function gi'^{x,Q'^) receives only gluonic contributions and is obtained as [4] 



g?'^{x,Q-') = 2el 



' 2tt 



ax y \y y 



(6) 



^The corresponding expressions in the case of additional weak boson exchange are given in [10]. 
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with 



p{3 - Az) - (1 - 2z) In 



1 + 



1-/3 



and AG(a;, Q^) the polarized gluon distribution. Here a denotes the threshold a 
eq the charge of the produced heavy quarks, and (5 is the cms velocity 



/3 = 



l-z 



(7) 

l + 4M2/g^ 
(8) 



of the final state quarks, with j3 e 
contribution to gi'^{x, Q"^) ^. Note that 

»l/a 



0, ^l-4M^x/[Q^{l-x)] ] . Eq. (6) defines the LO twist-2 



dzC^^^{z,M'Q,Q')=0 



(9) 



cf. [19], which leads to a positive and a negative branch of the structure function gi{x, Q^) in 
leading order for a positive definite polarized gluon density AG(a;,(5^), see Figure 2a. In this 
order also the first moment of g^^{x, Q^) vanishes, since the r.h.s. of (6) is a Mellin-convolution 
of two functions with the support of C^_^^ {z , Mq , Q"^) being z e [0, 1/a]. Eq. (9) also holds for 

the gluonic contribution to C'^^{z,Ml,Q'^) in the asymptotic limit Q"^ > M| in NLO [3]. For 
quarkonic initial states this relation does not hold, see [20]. 

The choice of the coUinear factorization leads to difficulties in deriving the correct twist-2 
terms in the case of structure functions which contain also twist-3 contributions in the limit 
of vanishing mass scales, as g2{x,Q'^) and g^ix^Q"^) in the case of electro-weak currents. This 
was extensively studied in the past [7, 8, 10] using the light-cone expansion and comparing the 
results to those being obtained in parton-model approaches'^. In Refs. [5-7] it was shown, that 
for fermionic contributions to the polarized structure functions the well-known results being 
obtained in the light-cone expansion, see e.g. [8], can be obtained if one refers to the covariant 
parton model, cf. [21]. This is due to the fact that the kinematic assumption (5) which neglects 
all parton momenta in the transverse direction is in confiict with the fact that the nucleon spin 
vector has both transverse and longitudinal components. This may have impact on predictions 
of relations between the moments of polarized structure functions, the Burkhardt-Cottingham 
sum rule [22-24] or integral relations between these functions. In the following we will therefore 
refer to a general orientation of the parton spin vector and use the covariant parton model instead 
of the coUinear approach. 



3 The general Case 

We now consider the polarization asymmetry of the hadronic tensor PF^^^ for a general nucleon 
spin vector 5"^ and general gluon spin vector s^, respectively. The nucleon spin vector obeys : 

S^P" = 0, 5^ = 5f + 5^ 
-Sf = (0;0,0,M) 

Si = M(0;cosx,sinx,0) . (10) 

^Note, that mass-scale effects may introduce twist 3 contributions to the structure function gi{x, Q^) as weU, 
as has been shown for target mass corrections in Ref. [10]. 

^For a review of ear her results see the comparison given in Ref. [8]. 
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The latter two relations hold in the nucleon rest frame. 

The polarization asymmetry of the hadronic tensor reads 



Wl^f{q,p, S) = j d'k [/+(p, k, S) - f_{p, k, S)] wl^\p, q, k, s) . 



The nucleon spin S is assumed to enter (11) linearly as usually the case in single photon- 
fermion interactions [25]. Here, k denotes the gluon 4-momentum. The gluon distribution 
functions f±{p,k,S) refer to opposite proton spin directions and are supposed to be twist-2 
parton distribution functions in the present paper. The sub-system hadronic tensor asymmetry 
w^/^\p, q, k, s) depends in addition on the virtual photon momentum exchanged, q, and the 



gluonic momentum- and spin vectors, obeys 

p.k 



^{p.kfk'^ - M^k^ 



:i2) 



with s.k = 0, k.k = —k"^. We consider general values of the gluon virtuality k'^, which is assumed 
to be sufficiently damped by the difference of the distribution functions f±{p, k, S) as k'^ — > oo, 
in order to keep (11) a well defined relation. Let us denote 



^J{p.kYk^ - M^k^ = N . 
The following relations are obtained : 



(13) 



S.s = —p.k S.k 

N 

^ ^ L k(i- ^^J^ 

q.k N \ p.k q.k 



We define now 



A/ = ^ (/.-/-) 



S.k 
M2 



f{p',p.k,e) 



(14) 



(15) 



and construct the sub-system hadronic tensor wj^\p, q, k, s) for single photon exchange and 
general values of the virtuality k^. The Lorcntz structure is determined by the Lcvi-Civita 
symbol contracted with two 4-vectors of the problem. wj^^\p, q, k, s) obeys the representation 



wii\p,q,k,s) 



le 



M p.k 
N q.k 



kf" 



k^ 



kY 

p.k 

k^ 



[9i+92]-q''k'' 



k'^p.q 
p.kq.k 



92 



p.k 



(16) 



The hadronic tensor is thus given by 
W(^\q,p,S) = ±-e,.^p j d^kJ{p\p.Kk')^ 



q'^k^ y 

k^p^k"^ 
p.k 



_ k'^p.q _ ^ 

9i + — 75-2 
q.kp.k , 



p.k 



idi + 92) 



(17) 
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Here gi = gi{q^, q.k, k!^) and v — v{q^, q.k, k^) denote the respective sub-system structure func- 
tions. The function v emerges for k"^ 0. 
For later use we rewrite (17) as 



W(^\q,p,S) = ^e,^^^ J d'k S.k U 



k^ 



dai dbi da2 db2 

+ p.q 



dp.k dq.k 



dp.k dq.k 



+ q'^^ 



das db^ ^ dui db^ 



dp.k dq.k dp.k dq.k 



dp.k dq.k J 



where 



dai 
dp.k 
daj 
dp.k 
dbi 
dq.k 
db2 
dq.k 
dh_ 
dq.k 
dbi 
dq.k 
dh_ 
dq.k 



f{p\p.k,k^) 
k' 



f{p -iP'k-, A; ), for i = 2 . . . 5 

'p.Ki 
q.k 



1 



{q.kf 



g2{q'^,q.k,k'^) 



9i{(f,(l-k,k'^) 
92{q^,q-k,k^) 

~'qk ^(^^'^-^'^^^ ■ 



(18) 

(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 



The functions Oj = ai{p'^,p.k, k'^) and bi = bi{q'^, q.k, k"^) will be used in Eq. (42). 

The tensor structure in Eq. (17) is yet different of that in (1). To determine the struc- 
ture functions 511,2(2;, Q^) a tensor decomposition in terms of the outer variables p,q and 5' is 
performed : 



W;,^\q,p, S) = -±^s,.^f,q''Sr \lt + it+P^Jl +P^Jl\ + -^^e^.apSrP^K-' , (26) 



with 



a r kf^k'^ ~ 



M2 



(27) 



-'2 



Jr, 



. k^k'^ k"^ p.q 

V 

{q.kyp.k 



-I 



f-g2^ A2g^^ + V + Cs/g" + D2 {p^q^ + (28) 



—I 



d^k 



d^k 



p.k q.k 
k^ k^ 
p.k q.k 

p.k q.k 



f-gi^E,p- + H,q- (29) 

f ■ g2 = E2P^ + H2q^ (30) 

f-v^ A^g''^ + V + Cw"g" + {p'^q^ + g V) • (31) 
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The contributions due to Bi^2, Ci,2, -^1,2, By and Dy vanish. Ay has to vanish because of current 
conservation. Therefore Wj^^ is given by 



Wj^^^ {q, p, S) = ■^e^^apq^S'' [Ai + A2] + -^e ^,^pqy^ S .q [L>i + L>2 + //i + //a + Cy] 



The coefficients ^1,2, -Di,2, -f^i,2 and read : 

J2(^ K\2 



^1 
A2 

D2 
Hi 

Cy 



I 



d^k 
d^k 



e ^ q\p.kf 



jp-k) 

2{q.k) ' 2{p.q)\q.k) (p.q) 

k^{p.q) ^ q^k'^ip.k) 



f -91 



2{q.kY{p.k) 2{p.q){q.kY (q.k) 
P 3q\p.kf 2{p.k) 

2{p.q){q.k) 2{p.q)'^{q.k) {p-qY 



f -92 



f -91 



k' 



3q^k\p.k) 



+ 



2k^ 



—I 



d^k 



2{p.k){q.k)'^ 2(p.qy{q.k)^ {p-q){q-k) 
f -91 



f -92 



d^k- 



{p.q){q.k) 
fc2 



{p.q){q.ky 
/■ 4 e{p.k) ^ 



We finally obtain the following representations for gi{x, Q^) and g2{x, Q"^) : 
gi{x,Q')+g2{x,Q') = J^i^i + M 



+ 



+ 



q (p.k) p.k 
2q.k ' 2q.k {p.qY p.q 

k^ P-Q _|_ P-k k"^ 



2 [q.k^p.k 2p.q[q.kY q.k 



d'^k 



2q.k p.q 



p.k 



f -91 
M2 



+ 



f -gi 

J ■ 92 

k^ 



^[Ai + A2] +^^[Di + D2 + H, + H2 + Cy] 



I 



^4i, P-(l 



q^ip.kY _|_ p.k 



q.k q.k {p.qY p.q 



q k p.k 



d^k 



[q.kY p.q 
q^ {p.kY 



q.k p.q 

Here the functions $i^2(^,P, ?) are finite as A;^ — > 0. 



+ p.k 



f ■92- 
f■9l 



p.k k"^ 
q.k p.q 



M2 



+ / d'^k 



f -91 



(32) 

(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 



(40) 



^2{Kp,q) . (41) 
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4 Representation of the structure functions in terms of a 
generating functional 

Equivalently to the representation of the polarized structure functions in the previous section one 
may represent them referring to a generating functional^. The hadronic tensor can be represented 
by functions of the form 



(42) 



where the functions Oj represent the parts depending on p and k and hi have a dependence on 
q and k only. As shown in the foregoing this separation is possible accounting for other factors, 
which do not depend on k. Here / belongs to the former and ^i, ^2, and v to the latter besides 
of scalar products. 

The partial derivatives of F by g and p can be represented by partial derivatives of the 
functions a and h as follows : 



dF 
dq'' 



dF 
dp"" 



dF dF 



dp.q 

ct^k a 
dF 



d^F 



dp.q 
j d^k b 
dF 



dq^ 

, dh ^ db 
dq.k dq'^ 
dF 



(43) 



dp.q 



da „ f^o, 

dp.k dp"^ 

d^F ^ d^F 
+ ^PaPxTT^ V2q„qx- 



(44) 



d'k 



k\k 



da db 



dq^dp.q 



d^F , d^F 

PaQX^, T^T+MaPx 



X^a- 



+ '2kxqa 



da db 



+ '2'Pxka 



' d{p.qy 
da db 



+ 4:Pxqa 



Qq2.Qp2 

da db 



dp.k dq.k dp.k dq^ dp"^ dq.k dp^ dq^ 

The following structures contribute to the polarized part of the hadronic tensor : 

da db 



.(45) 



^fiuaaQ ^x J d k k k - 
^^ivaaP S\ j d k k^k'^ 



dp.k dq.k 
da db 
dp.k dq.k 



^fMuaaP 



'gdF ^ . d^F 
S^^- h S.qp^ 



dp.q 
dF 

S"- + 2S.qq' 



dp.q 
-2q''Sx J 



d^k ¥ 



d{p.qY 

d{p.qy 
da db 
dp.k dq^ 



(46) 



■'livacr' 



d^kk^ 



da 
dp.k 



■ b = ef,^ap(fjPS.q 



dF 
dp.q ' 



(47) 
(48) 



which is given by 

Wl,t\<l.P.S) 



M2 



dFi 



+ p-q-. 



dFo 



d{p.q) ■ d{p.q) 



+ S.q p 



d^F, 
d{p.qy 



+ p.q 



d'F2 

d{p.qy 



*This case was studied for polarized light quarks in the presence of an on-shell condition in Ref. [5] before. 
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d{p.qY 
+ 2S.q 



d{p.q) dq^ 



S.k da^ dbs 
S.q d{p.q) dq^ 



The comparison with (1) yields the following expressions for the structure functions : 



gi{x,Q^)+g2{x,Q'') 



92{x,Q^) = 



M2 



p.q- 



dFi 

dp.q ' ^ ^ dp.q 



(49) 



(50) 



M2 



d^F, , d^F2 ^d\F^ + F,) , ^ a^Fs 

+ p-q— — ^ H — — — — \- 2- 



d{p.qy 



d{p.qy dip.qY 



dp.qdq^ 



4^ S.k da5 dh 
S.q dp.q dq^ 



d^k 



(51) 



We now separate the finite contributions to the polarized structure functions in the limit 
A;^ — > from those which vanish. We consider 



±{x[g,{x,Q') + g,{x,Q')]} = 



where x = —(f/{2p.q). Likewise one obtains 

d_ 

dx 



d^Fi d^F2 dF2 

+ p-q— — 7^ + 



d(p.qy 



d{p.qY d{p.q) 



(52) 



X— {(7i(x,g2)+^2(a:,Q')} = ^i(x,g2)+^2(a:,g') 

{p.qf " 



+ 



M2 



d{p.q)' 



+ p.q 



dFi 

d{p.qy ' d{p.q) 



+ 



(53) 



On the r.h.s. of (53) one may express the structure function g2{x.i Q^) inserting (51) which yields 



d 



-x—{g,{x,Q'')+g2{x,Q'')} = gi{x,Q^) - ^x,Q'') , 



with 



{P-Q? 
M2 



g(F3 + F4 - F2) 

d{p.q) 



+ 2 



d^F, 



d{p.q)dq^ 



2 f d^k ^'^ 
J S.q dp.k dq^ 



(54) 



(55) 



Let us investigate the structure of the function $(a;, Q^) more closely. To do this we refer to 
(46-48) from which follows 



dF, 
dp.q 



I 



,4, S.k dcii 
d k^^—r ■ b . 
S.q op.k 



One notices that 



^./, for k^2...^. 



(56) 



(57) 



dp. k p.k 

The remainder terms contributing to $(a;, Q^) result from (48). Using Eqs. (56,57) one thus 
concludes that $(a;, Q^) obeys the representation 



$(x,Q^) = / d'k (^^) cl>(p.k,q.ky,q',e) , 
where (p{p.k, q.k,p^, q^, k"^) is finite for k^ 0. 



(58) 
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5 The relation between gi{x, Q^) and g2{x, Q^) 



In (40, 41, 50, 51) the virtuality of the gluon field is revealed by power corrections in {k'^/Q'^y. 
These functions are not yet projected onto the twist-2 contributions.^ Any kind of partonic 
approach is only valid if the partonic virtualities obey 

\k^\ < . (59) 

This is an analogous condition to the requirement that the parton lifetime has to be much larger 
than the interaction time in the deeply inelastic scattering process for all associated infinite 
momentum frames, [27] . To project out the twist-2 contributions we refer to the coUinear basis, 
cf. [28], and expand these functions into a Taylor series in k'^ at k^ = zp^. Here p.p — and the 
lowest twist contribution is obtained in setting A;^ ^ in (40, 41, 50, 51). Note, however, that 
the corresponding expressions may contribute at higher twist due to the associated derivatives 
in k'^. We denote the twist-2 contributions to the structure functions gi{x, Q"^) by gY^{x, Q"^). 

Whereas the relation between the structure functions g^ and g^ is not easily seen from 
Eqs. (40, 41) ^, it can be directly obtained from (54), 



d_ 
dx 



^"(^,Q') + ^?(^,Q')1=^"(^,Q'), (60) 



the differential form of the Wandzura-Wilczek relation, see [9, 10]. The integral form can be 
obtained from (60) using the condition 

lim^f(x,g2)^0 (61) 

a;— >1 

which yields 

g'^{x, Q') = -g'^ix, Q') + C %'^{y, Q') . (62) 

Jx y 

At the level of twist-2 all structure functions depend on the same non-perturbative function 
/and sub-system structure function gi^\ Since the Wandzura-Wilczek relation holds as well 
for quarkonic initial states [6, 7] all higher order radiative corrections can be absorbed into g^^^^, 
respectively, with par = g, Qi denoting the quark species and the gluon, and holds thus in all 
orders. 

The validity of the Wandzura-Wilczek relation also in the case of gluonic operator matrix 
elements is in accordance with the observation of the general nature of this relation connecting 
vector operator matrix-elements with the associated scalar operator matrix-elements on the 
light cone, which was shown for fermionic fields in [13, 15]. 

To compare the heavy fiavor contributions to gi^2ix,Q^) to the usual parameterizations of 
these structure functions numerically we show the light flavor contributions to xgi 2{x,Q'^) in 
Figures la,b at leading order. For the parton distributions we refer to the parameterization [29]. 
Other recent parameterizations [30, 31] are well in accordance with [29] within the current 
experimental errors. The polarized gluon and sea-quark distributions have still a rather large 
uncertainty which can only be lowered by more precise data from future experiments. The present 
paranieterizatioiis \ycvc deri\Td assuming three light quarks through a fit to the data ignoring 



^Plain consideration of off-shell contributions in to hadronic structure functions may introduce unphysical 
contributions. One example are wrong target mass corrections as noted in [26] a long time ago. 
®For this the explicit dependence on {p.q) had to be known for Eq. (40). 
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any heavy flavor contribution. Part of the present gluon-distribution is thus corresponding to 
the contribution due to heavy flavors and one cannot add the distributions in Figures la, 2a or 
Figures lb, 2b in a simple way. 

As Figure 2a shows the effect of xg^'^lx, Q"^) is small at low scales but rises rapidly with 
and should be taken into account in future refined QCD analyses. Due to the larger quark 
mass and charge-suppression the 6-quark contribution to xgi{x, Q^) is smaller than that due to 
charm-quarks. While the distributions are positive for a; > 5 • 10~^ they change sign for lower 
values of X. In the present example the polarized gluon distribution is positive. The Wilson 
coefficient (7) changes sign, since 



1 



3 -In 



(63) 



\im C^^Q{z,Q') = /3>0, (64) 

and vanishes at threshold z = Q^/(l + 4Mq). Furthermore Eq. (9) holds. The oscillating 
behaviour of the Wilson coefficient implies lower relative heavy flavor contributions than in the 
unpolarized case. 

Figure lb shows xg2{x,Q'^) for the light flavors due to the Wandzura-Wilczek relation. 
xg2{x,Q^) is positive for small x values up to x ~ 10~^ and changes sign then. The inte- 
gral over the positive function xgi{x, Q^) is such larger than the subtraction term xgi{x. Q"^) in 
the former region while the subtraction term dominates in the latter region. The function has 
to have a positive and a negative branch since the Wandzura-Wilczek relation formally covers 
the Burkhardt-Cottingham relation'', 

dxg2ix, Q') = . (65) 

Jo 

Due to the additivity of twist-2 structure functions w.r.t. their parton contents the relation has 
to hold for each flavor separately. 

Similar conclusions can be drawn for xg2'^{x, Q^) shown in Figure 2b. xg2^{x, Q'^) is widely 
positive for smaller values of x < 2 • 10~^ and negative for larger x. Again we would like to stress 
that the present numerical results on xgi 2{x,Q^) are very sensitive on the polarized gluon 
distribution. 



6 Conclusions 

We calculated the heavy flavor contribution to g2{x, Q^) at leading order in the strong coupling 
constant using the covariant parton model for flnite values of the gluon virtuality k'^. The 
representation of the polarized structure functions gi2{x, Q"^) can be obtained applying a tensor 
decomposition. Furthermore, a generating functional in which the /c-dependent parton densities 
and coefficient functions are connected can be used to obtain a representation of the structure 
functions from which their possible relation can be derived. The twist-2 contributions are 
obtained in the limit k'^ <C Q'^. The functions gi^{x,Q^) obey the Wandzura-Wilczek relation 

^Note that the Wandzura-Wilczek relation is the analytic continuation from the positive moments, cf. e.g. 
[8, 10] in the local light cone expansion, where the 0th moment, which corresponds to the Burkhardt-Cottingham 
sum rule, does not contribute. 
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for a gluon induced process similar to earlier findings in fermion induced processes [5-10, 12, 13, IS- 
IS]. The absolute value of the relative numerical effect on both structure functions due to the 
LO heavy flavor contributions is about the same and may reach values of up to 5-10% in some 
kinematic ranges. ^ To make future QCD analyses of even more precise experimental data 
consistent at the level of the twist-2 contributions the heavy flavor distributions have to be 
taken into account. 
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Figure la: The light flavor contributions to the polarized structure function xg{{x,Q'^) in leading 
order (parameterization ISET=1 of [29], Aqcd = 203MeV) as a function of x and Q'^. Full line: 
g2 = 4GeV2; dashed line: ^ lOGeV^; dotted line: ^ loOGeV^; dash-dotted line: ^ 
lOOOGeV^ 
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Figure lb: The light flavor contributions to the polarized structure function xg2{x, Q'^). All conditions 
are the same as in Figure la. 
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Figure 2a: Heavy flavor contributions to the polarized structure function xgi{x^Q'^) in leading order 
(parameterization ISET=1 of [29], Aqcd — 203MeV) as a function of x and Q"^. Upper lines: 
gf{x,Q^) for rric = l.SGeV. Full line: = 4GeV^; dashed line: = lOGeV^; dotted line: = 
lOOGeV^; dash-dotted line: = lOOOGeV^ Lower lines: c/f(a;,Q2) for mt = 4.3GeV. Dotted 
line: = lOOGeV^; dash-dotted line: = lOOOGeV^. Figure 2b: Heavy flavor contributions 
to the polarized structure function xg2{x,Q^). 
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Figure 2b: Heavy flavor contributions to the polarized structure function xg2{x,Q'^). All conditions 
are the same as in Figure 2a. 
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